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We report on a simple scheme that may present a non-trivial geometric Zak phase {^zak) struc¬ 
ture, which is based on discrete-time quantum walk architectures. By detecting the Zak phase 
difference between two trajectories connecting adjacent Dirac points where the quasi-energy gap 
closes for opposite values of quasi-momentum (/c), it is possible to identify geometric invariants. 
These geometric invariants correspond to ~ ^ ~ ~ argue 

that this effect can be directly measured. 


1 .. INTRODUCTION 

Quantum physics attaches a phase to particles due to 
the complex nature of the Hilbert space. Phases arising 
during the quantum evolution of a particle can have dif¬ 
ferent origins. A type of geometric phase, the so-called 
Berry phase [T], can be ascribed to quantum particles 
which return adiabatically to their initial state, but re¬ 
member the path they took by storing this information 
on a geometric phase (ff>), defined as [Tl[2]: 

e*'*’ = (^/’inil^/’flnal)- (1-1) 

Geometric phases carry a number of implications: 
they modify material properties of solids, such as 
conductivity in graphene [3], they are responsible for 
the emergence of surface edge-states in topological in¬ 
sulators, whose surface electrons experience a geometric 
phase [4], they can modify the outcome of molecular 
chemical reactions [5], and could even have implications 
for quantum information technology, via the Majorana 
particle [6], or can bear close analogies to gauge field 
theories and differential geometry [7]. In this Letter, 
we present a simple system based on a discrete time 
quantum walk architecture, which presents a non-trival 
Berry phase structure on the torus, i.e., the Zak phase [8]. 

Discrete-time quantum walks (DTQWs) [9] offer a 
versatile platform for the exploration of a wide range 
of non-trivial geometric and topological phenomena 
(experiment) [TOHE], and (theory) p^UTS] . Further, 
QWs are roboust platforms for modelling a variety of dy¬ 
namical processes from excitation transfer in spin chains 
[Hilo] to energy transport in biological complexes m- 
They enable to study multi-path quantum inteference 
phenomena [22H25] . and can provide for a route to 
validation of quantum complexity [26l |27] , and universal 
quantum computing [28]. Moreover, multi-particle QWs 
warrant a powerfull tool for encoding information in an 
exponentially larger space, and for quantum simulations 
in biological, chemical and physical systems, in ID and 
2D geometries HHEgHal]. 

In this paper, we present a simple theoretical scheme 
for generation and detection of a non-trivial invariant ge¬ 


ometric phase structure in ID DTQW architectures. The 
basic step in the standard DTQW is given by a unitary 
evolution operator U{0) = TRfiiO), where RniO) is a ro¬ 
tation along an arbitrary direction n = (n^,, riy^riz)^ given 

by 

cos(^) — iuz sin(^) (m^, — Uy) sin(^) \ 
{irix + Uy) sin(^) cos(6>) + iUz sm{0) J ’ 

in the Pauli basis m- In this basis, the y-rotation is 
defined by a coin operator of the form mi- 

jy { cos(6») - sin(6») \ 

^ y sm{0) cos(^) J ' 

This is followed by a spin- or polarization-dependent 
translation T given by 

r = ^ |x + l)(:r| O \H){H\ + |x - l)(:r| ® \V){V\, 

X 

where H = (1,0)^ and V = (0,1)^. The evolution oper¬ 
ator for a discrete-time step is equivalent to that gener¬ 
ated by a Hamiltonian H{0), such that U{0) = 

= 1), with 

H{e) = f dk[E0{k)n{k).a] G \k){k\ 

J —TV 

and a the Pauli matrices, which readily reveals the spin- 
orbit coupling mechanism in the system. The quantum 
walk described by U{0) has been realized experimentally 
in a number of systems |32l [38ti4Q] . and has been 
shown to posses chiral symmetry, and display Dirac-like 
dispersion relation given by cos{Eo{k)) = cos{k) cos{0). 

Here, we present two different examples of non-trivial 
geometrical phase structure. The first DTQW proto¬ 
col consists of two consecutvie spin-dpendent transla¬ 
tions T and rotations R, such that the unitary step be¬ 
comes U{01,62) = TR{0 i)TR{ 02), as described in detail 
in m- The so-called “split-step” quantum walk, has 
been shown to possess a non-trivial topological landscape 
characterized by topological sectors with different topo¬ 
logical numbers, such as the winding number W = 0,1. 
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The dispersion relation for the split-step quantum walk 
results in m- 

cos{Eo^(l){k)) = cos{k) cos(^i) cos(^ 2 ) — sin(^i) sin(^ 2 )- 

The 3D-norm for decomposing the quantum walk 
Hamiltonian of the system in terms of Pauli matrices 
Hqw = E{k)n • a becomes [10]: 




sin(/c) sin( 0 i) cos(02) 
sin(£:6)^ ,02 (^)) 

cos(/c) sin( 0 i) cos(02)-|-sin(02) cos( 0 i) 

s\n{Ee^,e 2 (^)) 

— sin(A:) cos(02) cos( 0 i) 
s\n{E 0 ^^ 0 ^{k)) 


( 1 . 2 ) 


The dispersion relation and topological landscape for 
the split-step quantum walk was analyzed in detail in 
m- We now turn to our second example. 


The second example consists of two consecutive non¬ 
commuting rotations in the unitary step of the DTQW. 
The second rotation along the x-direct ion by an an¬ 
gle such that the unitarity step becomes 0) = 
TRx{(l))Ry{0)^ where Rx{(t>) is given in the same basis 
IHby: 


R,(^) = { "‘"W ) . 

Y 2 sm((/)) cos((/)) J 

The modified dispersion relation becomes: 

cos{Eo (f){k)) = cos{k) cos{0) cos(0) + sin(/c) sin(^) sin(0), 

(1.3) 

where we recover the Dirac-like dispersion relation for 
0 = 0, as expected. The 3D-norm for decomposing the 
Hamiltonian of the system in terms of Pauli matrices 
becomes: 

mX _ - cos(/c) sin(0) cos(g)+sin(/c) sin(^) cosjcj)) 

0,(f)\ ) sin(£;0,0(/c)) 

r^y (U\ _ cos(/c) sin(0) cos((/))+sin(/c) sin((/)) cos(0) 

^ 0 ,A'^> - siniEe,^ik)) 

/7 \ _ — sin(/c) cos(0) cos(0)+cos(/c) sin(0) sin(0) 

- sin(E 0 ,^(k)) • 

Dispersion relations for the DTQW with two consecu¬ 

tive non-commuting rotations within the first Brillouin 
zone are of the Dirac type (linear) at the set of gapless 
Dirac points, where the quasi-energy gap closes at 
E{k) = 0. The second rotation enables to close the gap 
at zero energy for complementary points, and allows 
to create a non-trivial geometric phase structure in the 
system. In particular, this system has a non-trivial 
phase diagram with a larger number of gapless points for 
different momenta as compared to the system consisting 
of a single rotation. We calculated analitically the 
gapless Dirac points for the system. Using basic trigono¬ 
metric considerations, it can be shown that the energy 
gap closes at 13 discrete points, for different values of 
quasi-momentum k. The phase diagram indicating the 
Dirac points where the gap closses for different momen¬ 
tum values is shown in Fig. 1. Squares correspond to 


|k|=Ti 



Ry-rotation (0) 

FIG. 1. Non-trivial phase diagram for the quantum walk 
with consecutive non-commuting rotations, indicating gap¬ 
less Dirac points where cuasi-energy gap closses for differ¬ 
ent values of quasi-momentum: Squares (k = 0), pentagons 
(|A;| = tt), romboids {k = -h7r/2), circles {k = —7r/2). 


Dirac points for /c = 0, circles correspond to Dirac points 
for k = —7r/2, romboids correspond to Dirac points for 
k = +7r/2, and pentagons correspond to Dirac points 
for \k\ = TT. This geometric structure in itself is novel 
and non-trivial and had not been studied in detail before. 


2.. ZAK PHASE CALCULATION 


We will now give expressions for the Zak Phase in two 
different scenarios. These scenarios are casted by the 
following hamiltonian 

H ^ UxCFx + riyay + (2.5) 

The hamiltonian to be described differ by a multiplying 
factor and by the expression of the But since the 
eigenvectors are the only quantities of interest for the 
present problem, the overall constants of this Hamilto¬ 
nian can be safely ignored. Now, our generic hamiltonian 
is given by the matrix 


^ nx-iuy \ 

\nx-\-iny -Uz J ' 

and has the following eigenvalues 

A = ±^Jnl + nl+nl 
The normalized eigenvectors then result 


>= 




UxEiriy 

2n2+2n2+2n2=p2n^ydi2+n2+n 

nzTy/n^+n^+n^ 

y2nJ+2n2_|_2n2=F2ras 



( 2 . 6 ) 


(2.7) 


( 2 . 8 ) 
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Note that the scaling rii Xrii does not affect the re¬ 
sult, as should be. This follows from the fact that two 
hamiltonians related by a constant have the same eigen¬ 
vectors. 

The overall Zak phase for the problem is to be consid¬ 
ered below is 




p7Tl2 

i / (< V^\dkV^ > + < V-\dkV- >)dk. (2.9) 

J — 'KI 2 


We will now apply these concepts to some specific exam¬ 
ples. 


2.1. Split-step Quantum Walk 

We first consider the split-step quantum walk. This 
corresponds to a quantum walk with unitary step give 
by U{ 0 i^ 02 ) = TR{ 0 i)TR{ 02 ), as proposed in [T^- In 
this example the normasl rii are of the following form 


”•@ 1.02 
0i,92 


{k) = 


y (k) 




sin(/c) sin( 0 i) cos(02) 

s\n{E 0 ^^e 2 (^)) 

cos(A:) sin( 0 i) cos(02)-|-sin(02) cos( 0 i) 

s\n{E 0 ^^e 2 (^)) 

— sin(/c) cos(02) cos( 0 i) 
s\n{E 0 ^^ 0 ^{k)) 


( 2 . 10 ) 


We consider the particular case that = 0. By taking 
one of the angle parameters such that Uz = 0 , it follows 
that the eigenvectors of the Hamiltonian are: 


1 / pi4>{k) 


tan 0 (/c) = —. ( 2 . 11 ) 


There are two choices for Uz = 0, which are = 0 or 
6>2 = 0. The phase is in both cases results in 


Z = Zx 


p7rl2 

^Z_=i dk<V+\dkV+> ( 2 . 12 ) 

»/—7r/2 


/*7r/2 

-hi/ dk <V-\dkV->= (I){-'k/2)-(I){'k/2), (2.13) 

J — TV 12 


f7r/2 
I-tv/2 

from where it follows that 


Z = 


tan( 6 > 2 ) 


tan(^i) 

A plot of the Zak phase is presented in Fig. 2 (a). 


(2.14) 


rriz = cos(/i:)c — sm{k)d (2.15) 

with 

a = sin(0) cos(6>) 

b = cos(0) sin(6>), (2.16) 

c = sin(0) sin(^), d = cos(0) cos(^). 

the angular functions defined above. The numerator Ni 
is given by 

Ni = Ux iriy = — exp(—iA:)(a — i6), (2-17) 

The remaining numerator N 2 is 

N 2 = ^z T -h = cos{k)c — sm{k)d 

^ a‘^ b‘^ cos‘^{k) + sin^(/c) — sm{2k)cd (2.18) 

On the other hand, the denominator D is reduced to 

D± = y2n| + 2n2 + 2nl T + ng + 

= (^a^ + 6^ + cos^(A:) + sin^(A:) — sin(2fc)cd 
^{cos{k)c — sm{k)d) 

1 

X \Ja? +h‘^ + cos'^{k) + d? sin^(fc) — sin(2A:)cdj . 


(2.19) 

By taking into account these expressions, one finds that 
the eigenvectors may be expressed simply as 


Ni 


V±>= Kf , <14|=(h^, ( 2 . 20 ) 


iV 2 ’ 
D± / 


El N 2 

Dl D± 


Then the calculation of the Zak phase 

^7r/2 
/-7r/2 


/>7r/2 

Z = Z^ + Z_=i dk< V^ldkYE > 

J—7rl2 


2.2. Quantum walk with non-commuting rotations 

The unitary step as described in the introduction re¬ 
sults in t/(d, 0) = TRx{(l))Ry{0). The norms rii are of the 
following form 

= — cos{k)a + sin(/c)5, riy = cos{k)b + sin(/c)a. 


p7Tl2 

+ i / dk<V-\dkV-> ( 2 . 21 ) 

—7r/2 


requires to know the following quantities 

r f m N 2 

Z±=ij 











4 


dkD± ]dk 


(|iViP + |A^2p) 
Di 


This expression can be simplified further. Due to (2.17 
it follows that the first term is real. However, an inspec 
tion of (2.18) shows that the last two terms are purel 


imaginary. Since the overall phase should be real, it fo] 
lows that these terms should cancel. This can be seen b 
taking into account that: 







D± = V\Ni\^ + \N2 


dk\Ni\^ = 0, (2.23, 



together with the fact that N 2 is real. Then 

2N2dkN2 


dkD± = 


2 Da 


(2.24) 


where (2.23) has been taken into account. Therefore 


•/(S' 


Z±=i ! [ + ^dkN2 


(|iViP + |A^2p) 


Di 


N2dkN2 ]dk, 


(2.25) 


but since D\ = |A^iP + |A( 2 p a simple calculation shows 
that the last two terms cancel each other. Thus the phase 
is 





NtdkN^ 

Dl 


dk. 


(2.26) 


By taking into account (2.19) the phases are expressed 
as 


Z± = 




{a^ + b^)dk 


Dl 


(2.27) 


We note that in this example the case = 0 is 
completely different than in the previous case, as it re¬ 
turns a trivial Zak phase Z = tt, since the k-dependence 
vanishes. We note that for this system the Zak phase 
landscape can be obtained by numerical integration. In 
particular, at the Diract points indicated in Figure 1, 
the Zak phase is not defined. 


A plot of the Zak phase ^zak is shown in Fig. 2, 
for parameter values 0i^2 = [~'7r,7r], and (j) = [—7r,7r]. 
(a) Zak phase for split-step quantum walk, given by 
the analytic expression Z = (b) Zak phase for 

quantum walk with non-commuting rotation obtained 
by numerical integration of expression Eq. 2.30. 


FIG. 2. (a) Non-trivial geometric Zak phase landscape for 
’’split-step” quantum walk, obtained analytically, (b) Non¬ 
trivial geometric Zak phase landscape for the quantum walk 
with non-commuting rotations, obtained by numeric integra¬ 
tion. 


the Brillouin zone. However, a geometric invariant 
can be defined in terms of the Zak phase difference, 
between two evolved states IV^finai)) which differ 

on a geometric overall phase only, starting from the 
same initial state (IV^ini)) at a saddle Zak point. In 
this case, the Zak phase difference can be written as 

('^finail'^finai) = Zak\. In [HI [43], experimental 

schemes for measuring the Zak phase difference for 
Bloch bands with two different dimerizations 
characterized by A/ci > 0 and A/c 2 < 0 where proposed, 
for [42] ultra-cold atoms, m photons in waveguides. 
In our system, we can define a geometric analogue of 
dimerization, given by two evolved trajectories (IV^finai)) 
with final parameter values associated with Dirac gapless 
points characterized by values of quasi-momentum k of 
opposite sign, resulting in a relative quasi-momentum 
differences Ak of opposite signs. These trajectories can 
be defined for all 13 Dirac points in the system since 
they are all analogue, up to arbitrary offsets which do 
not affect the calculation of the Zak phase difference. 
The two types of trajectories (IV^finai)) equivalent, up 
to an overall constant geometric phase. The appearance 
of a Zak phase difference of tt between the positive 
{Ak > 0) and negative {Ak < 0) branches, connecting 
singular points where the quasi-energy gap closes, 
signals a non-trivial geometric invariant structure in the 
system. On the other hand, for the case of trajectories 
connecting Dirac points between two positive {Ak > 0) 
or two n egative {Ak < 0) branches, the Zak phase 
difference should be zero. We note that in principle the 
dynamic phase does not cancel for each trajectory, but 
it does cancel for the difference between trajectories, as 
explained bellow. 


3.. GEOMETRIC INVARIANTS 

It is well known that the Zak phase is not a geometric 
invariant, since it depends on the choice of origin of 


The fact that the system aquires a different phase for the 
two trajectories (I'^finai)) corresponding to Dirac points 
for quasi-momenta of opposite sign can be readily seen 
by replacing the final parameter values associated with 
each Dirac and Zak saddle point in the DTQW unitary 
step. Let us consider the initial state I'lpini) and the start- 
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ing point for a trajectory at the center of phase diagram 
{0 = 0,0 = 0). We note that the exact initial state 
and trajectory is irrelevant, as long as the final states are 
equivalent up to an overall phase, since the Zak phase 
difference for our system only depends on the final states 
('^finail'^finai) ~ zak~^zak \. Wc cvolvc the systcm to 
either of the four adjacent Dirac points characterized by 
parameters {0 = 7r/2,0 = 7r/2), {0 = 7r/2,0 = —7r/2), 
{0 = —7r/2,0 = 7r/2), {0 = —7r/2,0 = —7r/2) by replac¬ 
ing these parameter values in the unitary step U{0,(l)) = 
TRx{(t))Ry{0). For 0 = (j) (positive branch) we find two 
Dirac points, which corresponds to A/c > 0, as indicated 
with romboid in Fig. 3. Evaluation of the product of ro¬ 
tations Rx{(t))Ry{0)^ for these parameter values, in a ma¬ 
trix which is proportional to cFz up to an overall phase 
g 27 r/ 2 Ar dependent on the step number N. This is the 
dynamical phase acquired by the system during its evo¬ 
lution. Therefore, up to an overall dynamical phase, the 
system returns to its initial state in the coin (or spin) sub¬ 
space, in addition to a translation in position to the site 
located at x = for initial state V(H) polarized 

(0(1), 1(0))^, due to the sequential translation operation 
. For the case of the trajectories evolved to adjacent 
Dirac points characterized by ^ = —cj) (negative branch), 
other two Dirac points are identified corresponding to 
A/c < 0, as indicated with circles in Fig. 3. The com- 



-Ti Ry-rotation (0) n 


FIG. 3. Geometric Zak phase difference \^zak ~ ^Zak\ 
trajectories connecting Dirac points with positive (negative) 
momentum difference Ak > (<)0 (positive(negative) branch), 
indicated by red (blue) dotted lines. The Zak phase difference 
is expected to be of the form ~ ^ between 

trajectories connecting points with opposite Ak and — 

~ ^ trajectories with Ak of the same sign (see 
text for details). Analogue geometric phase invariants can be 
defined for all other 12 singular points in the system where 
the quasi-energy gap closes. 


posed rotation for these parameters becomes 



This corresponds to exactly the same evolution with 
the same dynamical phase dependent on the number 
of steps but multiplied by an overall constant 

geometric phase Therefore, in this case the system 
is not only displaced by the translation operator to 
positions x = ±A, but in addition it aquires an overall 
constant geometric Zak phase of tt, as well as an overall 
dynamic phase which cancels out when measuring the 
phase difference Measurement of a Zak 

phase difference of either tt or 0 between trajectories 
(I'^finai)) corresponding to adjacent Dirac points would 
prove a non-trivial invariant geometric pattern in the 
system, which repeats itself over the entire Brillouin zone. 

Numerical simulations of the evolution of the probabil¬ 
ity distribution are shown in Figure 4, The overall phase 
(both dynamic and geometric) is of course not apparent 
in the probability distribution. The actual initial state 
and trajectory are of no concern, since the Zak phase 
difference only depends on the final states via 

(^finallV’final) = 

Fig. 4 (a) is the probability distribution for the 
initial state IV^ini) (step N = 0) located at x = 0, with 
righ-hand circular polarization l/\/2(l, i)^. Fig. 4 (b) is 
the probability distribution for an evolved state l^finai) 
with arbitrary parameters (6> = 7r/4, 0 = 0) after N = 7 
steps, which displays the typical ballistic delocalization 
associated with the Hadamard DTQW, Fig. 4(c) is the 
probability for the evolved state IV^gnai) adjacent 

Dirac point located at {0 = 7rl2,(l) = —7r/2) (negative 
branch) after N = 7 steps, where the system aquires an 
overall dynamic phase and a constant geometric 

Zak phase as explained above, and Fig. 4 (d) is the 
probability for the evolved state l^^nai) another adja¬ 
cent Dirac point located at {0 = 7r/2,0 = 7r/2) (positive 
branch) after N = 7 steps where the system acquires 
only an overall dynamical phase but no geometric 

Zak phase, according to our definition of the Brillouin 
zone. We note that the absolute Zak phase aquired in 
each trajectory is somewhat arbitrary and depends on 
the choice of origin, therefore it is not invariant. On 
the other hand, the phase difference picked up by the 
two trajectories described in Fig. 4 (c) and 4 (d) is 
geometry invariant and can be directly measured, this 
corresponds to a Zak phase difference \^zak~^Zak\ ~ 


4.. DISCUSSION 

A simple experimental scheme to measure the Zak phase 
difference can be envisioned. Starting from any Dirac 
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N=0 (0,4))=(n/4,O) N=7 



(0,4))=(n/2, -n/2) N=7 (0,4))=(ti/2, ti/2) N=7 



I4)fi„al'> 


0.5 

0.4 

(c)| 

1 1 

J 

1 4)“vn=Nn/2 1 

1 I 

1 

0.3 

ST 

0.2 

0.1 

0.0 


1 4)'>''"=Nn/2 1 

1 4)Z3k=0 1 


-10 -5 0 5 10 -10 -5 0 5 10 


FIG. 4. (a) Initial state {N = 0) located at a: = 0, with polar¬ 
ization l/\/2(l,z)^, (b) probability distribution after iV = 7 
steps with parameters (0 = 7r/4,0 = 0) corresponding to 
Hadamard qunatum walk, (c) probability distribution after 
iV = 7 steps with parameters {0 = 7r/2,0 = — 7r/2). With 
probability P = 1/2 the system is translated to the left at 
position X = —N and with probability P = 1/2 the system 
is translated to the right position x = N. During this evo¬ 
lution the system aquires a dynamical phase and a 

non-zero geometric Zak phase (d) probability distribu¬ 
tion after N = 7 steps for parameters (0 = '7r/2, 0 = '7r/2). 
With probability P = 1/2 the system is translated to the left 
at position x = —N and with probability P — Xjl the sys¬ 
tem is translated to the right at position x — N . During this 
evolution the system acquires a dynamical phase and 

zero geometric Zak phase. Note that the absolute Zak phase 
for each trajectory depends on the definition of the Brillouin 
zone. The Zak phase difference between these two trajectories 
~ = tt is a geometric invariant which can be 

measured. Equivalent geometric invariants can be measured 
for all the 13 Dirac points in the system. 


point in the parameter space, the system can be forced 
to evolve toward adjacent Dirac points through two 
different paths, corresponding to the positive (A/c > 0) 
branch and the negative (A/c < 0) branch. By shifting 
the origin of the Brillouin zone to the chosen Dirac 
point, the final parameter values for the negative path 
are characterized by [d = +(—)7r/2,0 = —(+)7r/2), 
while the final parameter values for the positive path 
correspond to {d = +(— )7r/2,0 = +(— )7r/2). For each 
path, a different geometric Zak phase should be aquired, 
such that = tt. This phase difference 

can be measured by recombining the two evolved states, 
in the case of photons by interferring the two evolved 
states via a Mach-Zehnder interferometer. On the other 
hand, for the case of two trajectories evolving along two 
positive or negative branches, the Zak phase difference 
should be zero = 0- ^ suitable scheme 

for detection of the Zak phase difference in a photonic 
system is described in [44]. 
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